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jj mE: 准确 估算 草地 净 初 级 生产 力 (Net Primary Productivity, NPP) 是 了 解 草地 生态 系统 碳 循环 过 程 和 科学 评估 
草地 对 气候 变化 响应 与 适应 机 制 的 关键 。 以 黄土 高 原 草地 生态 系统 为 研究 对 象 ,基于 1788 个 草地 生物 量 数据 及 19 


个 环境 因子 数据 (包含 
时 空 动态 ,并 估算 了 共享 社会 


气候 植被 .土壤 和 地 形 因子 ) ,利用 随机 森林 模型 模拟 了 2002 一 2020 年 黄土 高 原 草地 NPP 的 
经 济 路 径 (Shared Socioeconomic Pathway ,SSPs)4 个 未 来 气候 情景 下 ,黄土 高 原 草 地 
NPP 的 未 来 演变 趋势 。 结 果 表 明 :(1) 随机 森林 模型 的 模拟 精度 较 好 ,可 用 于 黄土 高 原 草 地 NPP 估 算 ;(2) 黄土 高 


在 陕西 关中 地 区 .甘肃 西部 地 区 以 及 山西 北部 地 区 


关键 词 : 草地 生 


草地 是 陆地 生态 系统 主要 组 成 部 分 之 一 , 约 占 
陆地 面积 30%0。 草 地 生态 系统 在 维持 生态 系统 功 
能 以 及 服务 社会 经 济 中 均 发 挥 重要 功能 ,包括 保护 
生物 多 样 性 .水土 保持 以 及 减缓 温室 效应 等 ”。 植 
被 净 初 级 生产 力 是 衡量 生态 系统 固 碳 能 力 的 重要 
指标 之 一 ”。 草 地 净 初 级 生产 力 (Net Primary Pro- 
duction,NPP) 不 仅 能 够 直接 反映 草地 生态 系统 生产 
力 的 高 低 , 而 且 是 评价 草地 生态 系统 矶 收文 及 调 市 
草地 生态 过 程 的 主要 因子 ”。 因 此 ,科学 评估 草地 
NPP 及 其 时 空 动态 对 草地 生态 系统 管理 及 草地 资源 
保护 至 关 重 要 。 

传统 野外 采样 估算 草地 NPP 的 方法 虽然 结 
较为 准确 ,但 由 于 其 需要 耗费 大 量 的 人 力 .物力 和 
财力 ,在 大 空间 尺度 草地 NPP 的 估算 中 较 少 使 用 。 
遥感 技术 和 生态 模型 的 发 展 为 此 提供 了 重要 的 方 


收 稿 日 期 : 2022-05-26; ”修订 日 期 : 2022-08-13 


原 草地 NPP 在 空间 上 整体 呈现 “东南 高 西北 低 ” 的 空间 分 布 特征 ,年 均值 为 276.55 g Com” ,其 中 陕西 关中 地 区 草地 
NPP 最 高 ;(3) 2002 一 2020 年 ,黄土 高 原 草地 NPP 总 体 呈 现 增加 趋势 ,其 中 55.01% 的 区 域 草地 NPP 增 加 ,主要 集中 
;(4) 在 气候 暖 湿 化 背景 下 
兽 加 趋势 ,其 中 SSP585 情 景 下 草地 NPP 增 加 最 多 ,SSP126 情 景 下 增加 最 少 。 利 用 随机 森林 模型 能 够 较 好 模拟 黄土 
高 原 草 地 NPP 时 空格 局 及 未 来 演变 趋势 ,为 黄土 高 原 草 地 生态 系统 保护 及 可 持续 发 展 提供 数据 支持 。 

态 系 统 ; 净 初 级 生产 力 ; 随机 森林 ; 共享 


,到 本 世纪 末 ,黄土 高 原 草 地 NPP 均 呈 


社会 经 济 路 径 ; 黄土 高 原 


法 和 技术 手段 。 遥 感 技术 通过 利用 遥感 指标 数据 
与 草地 生物 量 的 相关 关系 来 估算 草地 NPP ,如 马 文 
红 等 中 基于 NDVI 以 及 地 上 生物 量 数据 模拟 了 
1982—2006 年 中 国 北方 草地 生物 量 空 间 分 布 及 动 
态 ; 王 营 等 4 利用 EVI 与 草地 干 物质 产量 回归 模拟 
T 2006—2008 年 甘南 草地 NPP 变化 。 模 型 模拟 同 
样 是 大 空间 尺度 草地 NPP 估 算 的 常用 方法 ,如 Guo 
等 "利用 CASA 模型 模拟 了 2001 一 2018 年 内 蒙古 草 
Jii NPP; CE EAE Fl] FA CASA 模型 和 Miami 模 型 模 
拟 了 中 亚 草地 NPP;Zhang 45 * f] Hl BIOME-BGC 模 
型 模拟 了 呼伦贝尔 草原 NPP。 近 年 来 机 器 学 习 方法 
在 植被 生产 力 估算 中 亦 有 广泛 的 应 用 。Jung E 
利用 机 器 学 习 中 模型 树 集成 算法 模拟 了 全 球 总 初 
级 生产 力 ; Yao 等 由 利用 模型 树 集成 算法 模拟 了 
1982—2015 年 中 国 0.1° 空 间 分 辩 率 的 总 初级 生产 
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力 。 随 机 森林 (Random Forest,RF) 是 基于 决策 树 的 
集成 机 器 学 习 算 法 ,目前 已 被 应 用 于 全 球 呈 .我国 西 
北 地 区 "及 青藏 高 原 “ 不 同类 型 草地 NPP 的 估算 。 
黄土 高 原 是 世界 上 最 大 的 黄土 覆盖 区 ,气候 干 
F ,水土 流失 严重 ,生态 环境 极其 脆弱 ”“。 随 着 退 
耕 还 林 还 草 工 程 等 生态 工程 的 实施 ,黄土 高 原 植被 
变化 及 其 环境 效应 备 受 关注 ”。 目 前 ,在 黄土 高 原 
草地 NPP 时 空 动态 及 驱动 力 分 析 方 面 已 开展 了 一 
定 的 人 研究“”。 但 由 于 黄土 高 原 地 区 地 形 复杂 ,不 
同 研 究 结果 之 间 存 在 较 大 差异 。 在 全 球 变 暖 背 景 
下 ,黄土 高 原 草 地 NPP 未 来 的 演变 趋势 仍 不 清楚 。 
因此 ,本 研究 基于 气候 、 植 被 .土壤 和 地 形 因 子 , 利 
用 随机 森林 模型 模拟 了 2002 一 2020 年 黄土 高 原 草 
地 NPP 的 时 空 动态 ,并 基于 共享 社会 经 济 路 径 
(Shared Socioeconomic Pathway ,SSPs)4 个 未 来 气候 
情景 数据 , 佑 算 了 黄土 高 原 草 地 NPP 的 未 来 演变 趋 


势 , 以 期 为 黄土 高 原 草 地 生态 系统 管理 及 可 持续 发 
展 提供 科学 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

黄土 高 原 地 处 我 国 黄河 流域 中 部 ,地 理 位 置 位 
于 1009 54' ~114°33' E, 33?43' -41? 16' N, 总 面积 达 
63x10* km"? (图 1), 横 跨 陕 西 HEN TE LIV 青 
海 .河南 以 及 内 蒙古 7 个 省 份 ,是 我 国 四 大 高 原 之 
一 ,同时 也 是 世界 上 最 大 的 黄土 覆 善 区 。 黄 土 高 
原 包 括 了 干旱 SET REUS SE UTE URE ,年 降雨 
量 在 150-800 mm, ,海拔 在 200~3000 m ,不 同 区 域 在 
气候 上 有 着 明显 差异 2 。 黄 土 高 原 的 植被 类 型 以 
森林 .草地 以 及 灌木 地 为 主 , 草 地 面积 约 占 总 面积 
60957, 

1.2 数据 获取 及 处 理 

通过 文献 .国家 地 球 系统 科学 数据 中 心 数据 以 
及 GEDI 生 物 量 数据 (https://daac.ornl.gov/) 获 取 黄 土 
高 原 草 地 地 上 生物 量 数据 1788 个 ,根据 根 冠 比 估 算 
地 下 生物 量 , 使 用 草地 生物 量 乘 以 转换 系数 (0.475) 
估算 草地 NPP(g Cm”), 作 为 实测 NPP 数据 。 

为 训练 模型 ,选取 19 个 特征 参数 包括 土壤 理化 
因子 (5) 、 地 形 因 子 (7) .气象 因子 (4) 和 植被 因子 
(8B)4 类 。 气 象 因 子 中 温度 和 降水 数据 来 源 于 国家 
地 球 系 统 科 学 数据 中 心 - 黄 土 高 原 分 中 心 (http:// 
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图 1 研究 区 概况 及 采样 点 分 布 
Fig. 1 Study area and sampling points 


loess.geodata.en ) ,根据 月 平均 温度 数据 合成 年 平均 
温度 (TEM) 数 据 以 及 年 生长 季 平 均 温度 (TEM410) 数 
据 , 根 据 月 最 低温 度数 据 合成 为 年 最 低温 度 (TMN) 
数据 ,根据 逐 月 降水 量 数据 合成 年 降水 量 (PRE) 数 
据 以 及 年 生长 季 降 水 量 (PRE,1) 数 据 中 。 数 字 高 程 
模型 (DEM) 数 据 来 自 NASA ÉJ SRTM 数据 (http:// 
srtm.csi.cgiar.org/) ,坡度 (Slope) 数 据 由 DEM 数据 通 
TE ArcGIS 计算 得 来 。 土 壤 理化 因子 数据 均 来 自 世 
界 土壤 数据 库 (https://data.isric.org/)。 日 光 诱 导 叶 
绿 素 荧光 (SIF) 数 据 来 自 GOSIFv2 数据 "(https:// 
globalecology.unh.edu ) ,植被 有 效 光 合 辐 射 吸收 比例 
(FPAR ) 数 据 、 归 一 化 植被 指数 (NDVI) 以 及 燕 散 量 
(ET) 来 自 MODIS 数据 产品 (https:/modis.gsfc.nasa. 
未 来 降雨 量 和 温度 数据 来 源 于 WorldClim 
(https://www.worldclim.org/)。 黄 土 高 原 空间 范围 数 
据 以 及 黄土 高 原 草 地 分 布 数据 来 自 中 国 科学 院 资 
源 环境 与 数据 中 心 (https:Wwww.resdc.cn/) 空 间 分 辩 
率 为 1km。 以 上 数据 分 辩 率 均 通 过 ArcGIS FE 
为 1 km, 采 用 D_WGS_1984 投 影 。 本 文 所 有 环境 因 
子 见 表 1。 

SSPs 情景 在 典型 浓度 路 径 (RCPs) 情 景 基 础 上 
发 展 而 来 ,用 于 定量 描述 气候 变化 与 社会 经 济 发 展 
路 径 之 间 的 关系 ,反映 未 来 社会 面临 的 气候 变化 适 
应 和 减缓 挑战 。 本 研究 分 别 选择 SSP126、 
SSP245, SSP370 和 SSP585 4 个 情景 下 ,ACCESS- 
ESM1-5 ,BCC-CSM2-MR .CanESMS , CMCC-ESM? 和 
CNRM-CM6-1, 共 计 5 个 气候 系统 模式 2021—2040 
年 (2030s) , 2041—2060 4E (2050s) , 2061—2080 年 
(2070s) ,2081—2100 年 (2090s)4 个 时 间 段 的 降水 


gov/) 。 
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R1 随机 森林 CRF) 中 使 用 的 环境 因子 


Tab.1 Environmental factors in RF 


ESI 环境 因子 解释 分 辩 率 数据 来 源 
CLAY S 0-30 em SUL E 250m https://data.isric.org/ 
S CLAY T 30-100 cm ffr rtt 250m 
S SILT S 0-30 em 粉 粒 含量 250 m 
S SILT T 30-100 em 粉 粒 含 量 250 m 
S SAND S 0-30 em 砂粒 含量 250 m 
S SAND T 30-100 cm 砂粒 含量 250 m 
S SOC S 0-30 cm icd po 250 m 
S SOC T 30-100 cm 土壤 有 机 碳 含 250m 
T DEM 数字 高 程 模型 30 m https://srtm.csi.cgiar.org/ 
T Slope 坡度 30m 
A TEM 年 均 温 (2002 一 2020) 1000 m http://loess.geodata.cn 
A TEM 4—10 Fili (2002—2020) 1000 m 
A TMN 年 最 低温 度 (2002 一 2020) 1000 m 
A PRE 年 降水 量 (2002 一 2020 ) 1000 m 
4 PREsi 4—10 H 7k Œ (2002—2020) 1000 m 
A ET Zi wt (2002—2020) 500 m https://modis.gsfc.nasa.gov 
B SIF 日 光 诱 导 叶 绿 素 痰 光 (2002 一 2020) 0.05° hitps://globalecology.unh.edu 
B FPAR 植被 有 效 光合 辐射 吸收 比例 (2002 一 2020 ) 500 m https://modis.gsfc.nasa.gov 
B NDVI 归 一 化 植被 指数 (2002 一 2020 ) 1000 m 
注 :$ 表 示 土 壤 理 化 因子 ;7 表示 地 形 因 子 ;4 表示 气象 因子 ;有 表示 植被 因子 。 


量 和 温度 作为 输入 数据 ,利用 RF 模 型 预测 未 来 不 
同 气 候 情 景 下 黄土 高 原 草 地 NPP 的 演变 趋势 。RF 
模型 中 的 其 他 环境 因子 仍 使 用 2020 年 数据 。 

1.3 研究 方法 

1.3.1 建 模 方法 本 人 研究 采用 RF 进行 建 模 ,RF 是 
基于 分 类 和 回归 树 (Classification and Regression 
Tree,CART) 而 开发 的 集成 算法 ,可 用 于 分 类 问题 以 
及 回归 问题 了””。RF 的 主要 原理 为 通过 结合 多 个 
决策 树 ,并 平均 其 结果 使 得 决策 树 泛 化 误差 收敛 从 
而 产生 更 好 的 预测 结果 "小 。RF 可 广泛 用 于 预测 问 
题 , 并 且 应 用 简单 ,无 需 复 杂 地 调整 参数 ,通过 调整 
RF 中 回归 树 的 数量 以 及 决策 树 每 处 节点 预测 变量 
的 数量 来 优化 模型 ,提高 预测 精度 ”。 此 外 ,RF 具 
有 处 理 高 维特 征 空间 数据 的 能 力 ,在 遥感 应 用 中 已 
有 较 多 研究 ”。RF 模 型 性 能 通过 决策 系数 ( 尼 ) 以 
及 均 方 根 误差 (RMSE ) 评 价 其 拟 合 能 力 以 及 预测 性 
能 。 本 研究 使 用 R 4.0.5 进行 随机 森林 建 模 .参数 优 
化 .预测 以 及 空间 可 视 化 。 

1.3.2 Theil-Sen Median 趋势 分 析 和 Mann-Kendall 检 
XE Theil-Sen Median 趋势 分 析 是 一 种 稳健 的 非 参 


数 统计 趋势 计算 方法 3 ,计算 公式 如 下 : 


HN 
Sen = Median j 


slope 


=y 
= ) i,j=1,2, +, N (1) 


式 中 : x, Al x, 为 样本 数据 值 。 当 Sen gp > 0 时 ,表明 
该 数据 在 此 时 间 序 列 呈 增长 趋势 ,反之 则 表示 呈 下 
降 趋势 。 

Mann-Kendall 检 验方 法 是 一 种 非 参 数 统 计 检 验 
方法 ,用 于 判断 趋势 的 显著 性 ”。 计 算 公 式 如 下 : 


S= > > sign(x, - x,) (2) 


i=1j=itl 
Sel. 350 
var(s) 
A 5x 0; S=0 (3) 
S+1 g59 
var(s) 
P=21 - cdt(| zs |) (4) 


Cp : S X Mann- Kendall 检验 需要 构建 的 统计 量 ， 
Zur 是 标准 正 态 变量 。 本 文 在 置信 度 水 平 c=0.05 上 
来 检验 NPP 变化 趋势 的 显著 性 。 当 P 小 于 0.05 时 表 
示 NPP 变 化 显著 ,反之 则 表示 不 显著 。 
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120 干旱 区 研究 4035 
进行 验证 。 结 果 表 明 ,RF 的 训练 集 R?,, 20.97 , RM- 
2 结果 与 分 析 SE 4727.34 g Csm2; 进 行 十 折 交 叉 验证 后 的 测试 集 
R? =0.84, RMSE..=59.73 g Cnm-2( 图 3)。 模 型 精度 
21 NPP 与 环境 因子 相关 性 分 析 d. KA MS E 
l PES 较 高 ,可 以 用 于 黄土 高 原 地 区 草地 NPP 的 模拟 。 
在 植被 因子 中 ,NPP 与 SIF 相关 性 最 强 (r= 23 草地 NPP 空 间 动态 
0.82) ,其 次 为 NDVI 和 FPAR (图 2a)。 在 气象 因子 ~ ee ak NE gie eps i RE 
中 ,NPP 与 ET 相关 性 最 高 (r=0.91), 与 PREswo、PRE、 。 质 性 (图 和，2002_2020 年 黄土 高 原 草 地 年 平均 
TMN 相关 性 也 较 强 ,而 与 TEM 以 及 TEM。w 相 关 性 。 NPp 为 276.55g Com”. 黄土 高 原 草地 NPP 的 分 布 
65359, NPP 与 土壤 因子 整体 相关 性 较 弱 ,其 中 NPP 呈现 从 南 向 北 逐 渐 递减 的 趋势 .其 中 黄土 高 原 西北 
与 SAND_S 相关 性 最 强 (r2-0.27), H SOC T 相关 性 地 区 草地 NPP 最 低 , 主要 位 于 陕西 北部 宁夏 以 及 
最 弱 (r=0.11)。 整 体 而 言 ,NPP 与 气象 因子 .植被 因 甘肃 兰州 地 区 附近 ,最 低 值 为 76.52 g C «m^; ERR 
子 相 关 性 较 强 ,而 与 气象 因子 中 的 TEM All TEM uo, NPP 最 高 的 区 域 主要 分 布 在 黄土 高 原 南 部 地 区 ,如 
土壤 因子 以 及 地 形 因子 相关 性 弱 。 而 在 随机 森林 ”陕西 关中 地 区 ,最 高 值 为 634.86 g Com; TEE 
模型 因子 重要 性 中 (图 2b) ,对 模型 估算 NPP 影 响 最 高 原 东 部 地 区 草地 分 布 较 少 ,草地 NPP 约 为 450 g 
大 的 是 ET。PRE .PRE 以 及 植被 因子 中 的 SF 对 Cem 
NPP 估 算 影 响 也 较 大 。 结 合 相 关 性 分 析 和 重要 性 排 2002 一 2020 年 ,黄土 高 原 草 地 NPP 呈现 总 体 增 
FF. ,影响 黄土 高 原 草 地 NPP 的 环境 因子 主要 为 植被 加 的 趋势 (图 5), 其 中 16.59% 的 地 区 明显 著 增加 趋 
因子 以 及 PRE PRE ww 和 ET。 势 ,38.52% 的 地 区 旦 轻微 增加 趋势 ,主要 分 布 在 黄 
2.2 RF 模型 精度 评价 土 高 原 东部 以 及 南部 ;1.89% 的 地 区 草地 NPP 变化 
通过 多 种 植被 气候、 土壤 及 地 形 因子 建立 随 ” 不 显著 ;35.57% 的 地 区 明 轻 微 减少 趋势,7.43% 的 地 
机 森林 模型 ,利用 实测 数据 对 预测 的 草地 NPP 结 区 呈 显 著 减 少 趋势 ,主要 分 布 在 黄土 高 原 西 北 地 区 。 
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注 :NPP 为 净 初 级 生产 力 ;SI 为 日 光 诱 导 叶 绿 素 荧 光 ;FPAR 为 植被 有 效 光合 辐射 吸收 比例 ;NDVI 为 归 一 化 植被 指数 ;ET 为 蒸 散 量 ;PRE 为 年 
降水 量 ;PRE, 为 4 一 10 月 降水 量 ;TMN 为 年 最 低温 度 ;TEM 为 年 均 温 ;TEM4w 为 4 一 10 月 年 均 温 ;DEM 为 数字 高 程 模型 ;Slope 为 坡度 ;CLAY_S 
240-30 cm 黏 粒 含 量 ;CLAY _T 为 30~100 em RLE Ht ;SAND_S H 0-30 cm 砂粒 含量 ;SAND_T 为 30~100 cm 砂粒 含量 ;SILT_S 为 0~30 em BPEL 


含量 ;SILT T 2930-100 em 粉 粒 含 量 ;SOC_S 为 0~30 cm J 


图 2 NPP SHAS 


FE: 壤 有 机 碳 含量 ;SOC_T 为 30~100 em 土壤 有 机 碳 含 量 。 
E(a) 及 随机 森林 因子 重要 性 (b) 


Fig.2 The correlations between NPP and environmental factors (a), and the importance of random forest factors (b) 
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um a 2.4 未 来 不 同 气候 情景 下 黄土 高 原 地 区 草地 NPP 
。RF 测 试 集 Et 演变 趋势 
~ 50 到 本 世纪 末 ,4 个 未 来 气候 情景 下 黄土 高 原 草 
5 - 地 NPP 均 呈现 增加 的 趋势 ,其 中 SSP126 情景 草地 
= 400 NPP 增 加 幅度 最 少 ,为 120 g Cm? (103); SSP585 
z per 情景 增幅 最 大 ,为 8.10 g C-m?* (10a), SSP245 Ail 
Tomo grey SSP370 情 景 ,草地 NPP 的 增 速 分 别 为 3.90 g Com? 
E (10a) 814.00 g Cem? (102) (Fl 6), 2E HE NPP J 
olkoon i = 人 一， 加 的 区 域 主要 集中 在 黄土 高 原 南 部 地 区 ,SSP585 中 
“ ee 草地 NPP 在 黄土 高 原 南部 地 区 增加 最 多 。 而 草地 
NPP 下 降 的 区 域 主要 集中 在 黄土 高 原 西部 地 区 ,其 
图 3 随机 森林 模型 精度 验证 中 在 SSP126 中 草地 NPP 将 下 降 最 多 。 
Fig.3 Accuracy test of Random Forest model 
3 讨论 
100°E 105°E 110°E 115°E 
l | 3.1 随机 森林 模型 估算 草地 NPP 适 用 性 
相 比 其 他 NPP 佑 算 模型 ,随机 森林 的 优势 在 于 
£ J£ 可 以 更 准确 地 描述 NPP 与 其 他 环境 因子 之 间 的 非 线 
` ` 性 复杂 关系 , 且 随 机 森林 有 较 好 抗 过 拟 合 能 力 *™ 
本 研究 利用 idi 导 到 黄土 高 原 草 地 的 平 
2 z FJ NPP(276.55 g C.m2?)。 刘 洋洋 等 基于 CASA 模型 
a ` 1% 模拟 的 2000—2015 one. NPP Jj 
uiri De MM cm 202.93 g Com, IA RT AE 55 f JE RC A f e E] 
Dod rac M Bc 有 关 "”。 魏 建 洲 基 于 MOD17A3 NPP 数 据 模拟 的 黄 
"E is cai moe 土 高 原 草地 NPP 在 208.04~262.66 g C -m° zt li] ,其 
图 4 2002 一 2020 年 黄土 高 原 草地 NPP 空间 分 布 高 值 与 本 研究 结果 较为 接近 中 。Li 等 利用 CASA 模 
Fig. 4 Spatial pattern of grassland NPP in the Loess Plateau 型 分 别 模拟 了 石 羊 河 XAR 国 西北 干旱 半 干 早 
i E 区 草地 年 均 NPP, 其 结果 分 别 为 236.99 g C- mM 
— - - nee 277.97 g Cm sr” 。 相 比 于 其 他 研究 结果 ,本 研究 
TUR 模拟 草地 NPP 结 果 偏 大 ,一 方面 在 于 不 同 研 究 中 所 
国 显著 增加 使 用 的 草地 分 布 数据 不 同 ,如 刘洋 洋 使 用 的 草地 分 
al m | z 布 数据 来 自 Global Land Cover 2000 数据 产品 "”, 本 
| nn 3 研究 所 使 用 草地 分 布 数据 来 源 于 中 科 院 土地 利用 
数据 产品 ; 另 一 方面 基于 随机 森林 模型 模拟 草地 
NPP, 受 训练 数据 影响 很 大 ,本 研究 用 于 训练 模型 的 
E E 实测 草地 NPP 数 据 主要 集中 在 黄土 高 原 西部 ,东部 
m 地 区 的 样本 点 数量 相对 较 少 。 但 总 体 而 言 ,本 研究 
| l eee 利用 随机 森林 模拟 的 黄土 高 原 草 地 NPP 结 果 较 为 
100°E 105°E 110°E 115°E 可 信 。 
图 5 2002 一 2020 年 黄土 高 原 草地 NPP 空 间 变 化 32 黄土 高 原 草地 NPP 的 动态 变化 
Fig. 5 Spatial dynamic of grassland NPP in the Loess AME ZIM 2002—2020 E E + A Jg c Hb, NPP 5 


Plateau during 2002-2020 ISAIA, 3x 5j t A XE A DC BUE 2 SH FE gi 


202303.001 12v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


128 干旱 区 研究 40 卷 
100°E 104°E 108°E 112°E 100°E 104°E 108°E 112°E 
d d z 
e (a) SSP126 时 ien (b) SSP245 S 
El e g Z 
E & m a 
至 上 0 ZZ Z 
$45 & d a 
100°E 104°E 108°E 112°E 100°E 104°E 108°E 112°E 
100°E 104°E 108°E 112°E 100°E 104°E 108°E 112°E 
Z r T T T Z j T T T Z, 
Q | (c) SSP370 $ Q T (dSSP585 $ 
名 Z £ 名 
ea cea 加 
Z ZZ Z 
E! SA a 
100°E 104°E 108°E 112°E 100°E 104°E 108°E 112°E 


NPP/(g C-m?:2?) 


EEEESEERARBSO—] — 1 1 1 — —1 a E 
100 150 200 250 300 350 400 450 500 550 600 


图 6 不 同 气候 情景 下 2090s 黄 土 高 原 草地 NPP 空 间 分 布 


Fig.6 Spatial pattern of grassland NPP in the Loess Plateau in 2090s under different climate scenarios 


由 于 黄土 高 原 气 候 干 旱 , 降 水 是 主要 的 植被 生长 限 
制 因 子 " 外 。2002 一 2020 年 黄土 高 原 降 水 量 呈 增加 
趋势 (图 7) ,促进 了 黄土 高 原 草 地 NPP 的 增加。 刘 
洋洋 等 ”也 认为 黄土 高 原 降 水 增加 导致 的 湿润 环 
境 为 草地 NPP 增 加 提供 了 有 利 条 件 。Gang 等 ”的 
研究 结果 同样 表明 气候 变化 对 黄土 高 原 草地 NPP 
增加 有 着 促进 作用 。 从 空间 变化 上 来 看 ,黄土 高 原 
西部 和 南部 草地 NPP 增 加 ,与 刘铮 等 3 的 研究 结 
相似 。 这 些 地 区 位 于 半 干 旱 以 及 半 湿 润 地 区 ,降水 
量 的 增加 导致 了 NPP 的 快速 积累 。 此 外 , 自 1999 
年 以 来 ,退耕 还 林 还 草 工程 通过 围栏 封 育 、 封 山 禁 
牧 等 措施 使 草地 生态 环境 得 到 了 较 大 改善 ,是 导致 
黄土 高 原 草地 NPP 增 加 的 重要 因素 “2 。 到 本 世纪 
末 ,4 个 情景 中 黄土 高 原 草 地 NPP 仍 呈 增 加 趋势 , 原 
因 在 于 未 来 黄土 高 原 地 区 降水 量 和 最 低温 度 将 继 
续 增加 , 暧 湿 化 的 气候 条 件 有 利于 草地 NPP 的 积 
累 ,尤其 是 在 SSP585 情 景 中 。 
3.3 RF 估算 结果 的 不 确定 性 分 析 

本 人 研究 使 用 RF 模型 模拟 黄土 高 原 草地 NPP 精 
度 较 高 ,但 仍 存 在 一 定 的 不 确定 性 。 首 先 , RF 模型 
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7 2002 一 2020 年 黄土 高 原 年 降水 量变 化 趋势 
Fig.7 Temporal trend of mean annual precipitation in the 
Loess Plateau from 2002 to 2020 


性 能 受 数据 集 数量 及 质量 影响 较 大 ” ,黄土 高 原 区 
域 地 形 复 杂 ,本 文 使 用 1 km 的 空间 分 辩 率 数据 模拟 
草地 NPP 导致 部 分 区 域 草 地 NPP 模 拟 结果 精度 较 
低 。 其 次 ,RF 模型 作为 经 验 模型 ,受训 练 数据 质量 
以 及 数量 影响 较 大 ,本 人 研究 草地 NPP 实 测 数 据 空间 
分 布 不 均匀 ,在 黄土 高 原 西 北 地 区 以 及 西南 地 区 较 
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多 ,而 其 余地 方 较 少 ,未 来 应 考虑 使 用 更 多 的 实测 
数据 进行 模拟 。 此 外 ,RF 模型 属于 黑箱 模型 ,通过 
oe 定 合 理 模 拟 结果 ,但 无 法 控制 模型 内 


3 BAT 


部 运行 


4 结论 


nj 


本 研究 基于 野外 实测 数据 ,以 植被 ` 气 候 EE 
和 地 形 因子 作为 环境 因子 ,利用 随机 森林 模型 模拟 
了 黄土 高 原 草 地 NPP 的 时 空 动态 ,并 估算 了 未 来 不 
同 气 候 情 景 下 ,草地 NPP 的 演变 趋势 ,得 出 了 以 下 
主要 结论 

(1) 2002 一 2020 年 ,黄土 高 原 草 地 平均 NPP 为 
276.55 g C.m2.a, 空 间 上 呈 ”“ 东 南 高 西北 低 ” 的 分 布 
特征 。 近 19 a 黄 土 高 原 草 地 NPP 整 体 呈 现 上 升 趋 
势 ,主要 得 益 于 20 世 纪 以 来 黄土 高 原 地 区 降雨 量 的 
增加 及 “退耕 还 林 还 草 ” 等 一 系列 生态 工程 。 

(2) 到 本 世纪 末 ,在 上 暖 湿 化 的 气候 背景 下 ,黄土 
高 原 草 地 NPP 均 将 继续 保持 增加 的 趋势 ,其 中 在 
SSP585 情景 中 将 增加 最 多 。 

(3) 基于 气候 环境 .土壤 和 地 形 因 子 构建 了 黄 
土 高 原 草地 NPP 随机 森林 模型 ,尚未 包含 人 为 活动 
指标 因子 ,构建 综合 气候 环境 .土壤 和 人 为 因素 等 
多 因子 的 要 素 体系 ,将 会 进一步 提高 草地 植被 生产 
力 的 估算 精度 。 
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Simulation of spatial pattern and future trends of grassland net primary 
productivity in the Loess Plateau based on random forest model 
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Abstract: The accurate estimation of grassland net primary productivity (NPP) is crucial to understanding the 
carbon cycle of grassland ecosystems and their adaption to climate change. Based on 1788 grassland biomass 
data and 19 environmental factors (climate, vegetation, soil, and topographic factors), we simulated the 
spatiotemporal dynamics of grassland NPP in the Loess Plateau from 2002 to 2020 using the random forest (RF) 
model. The future trends of grassland NPP under four future climate scenarios of shared socioeconomic pathways 
were estimated. Results showed that (1) the RF model had a good accuracy, which indicated that RF can be used 
to estimate grassland NPP in the Loess Plateau; (2) grassland NPP in the Loess Plateau exhibited a “high in 
southeastern and low in northwestern” pattern, with a mean value of 276.55 g C* m ^-a '. The highest grassland 
NPP was observed in Guanzhong Plain of Shaanxi; (3) the grassland NPP in the Loess Plateau showed an overall 
increasing trend during 2002-2020. Regions experiencing an increase in grassland NPP accounted for 55.01% of 
the total land area, which is mainly located in Guanzhong Plain, western Gansu, and northern Shanxi; (4) under 
the wetter and warmer climate, grassland NPP in the Loess Plateau will continually increase by the end of this 
century. Grassland NPP will increase the most under the SSP585 scenario and the least under the SSP126 
scenario. RF can be used to simulate the temporal and spatial trends of grassland NPP in the Loess Plateau. The 
results provide data support for the protection and sustainable development of grassland ecosystem in the Loess 
Plateau. 


Keywords: grassland ecosystems; net primary productivity (NPP); random forest; Shared Socioeconomic Path- 
way (SSPs); Loess Plateau 


